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Abstract 

The temporo-spatial relationship between the three organs (lung, spleen and lymph node) involved during the initial stages 
of Mycobacterium tuberculosis infection has been poorly studied. As such, we performed an experimental study to evaluate 
the bacillary load in each organ after aerosol or intravenous infection and developed a mathematical approach using the 
data obtained in order to extract conclusions. The results showed that higher bacillary doses result in an earlier IFN-y 
response, that a certain bacillary load (BL) needs to be reached to trigger the IFN-y response, and that control of the BL is 
not immediate after onset of the IFN-y response, which might be a consequence of the spatial dimension. This study may 
have an important impact when it comes to designing new vaccine candidates as it suggests that triggering an earlier IFN-y 
response might not guarantee good infection control, and therefore that additional properties should be considered for 
these candidates. 
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Introduction 

The ability to subvert the host's immune system by blocking its 
ability to detect the pathogen is considered to be one of the main 
characteristics of Mycobacterium tuberculosis infection [1]. Different 
authors have elaborated several hypotheses to explain this 
phenomenon, most of which focus on the ability of M. tuberculosis 
to avoid apoptosis [2] and block antigen presentation [3,4] or 
lysosomal degradation [5], amongst others. However, these 
hypotheses are based on mechanisms studied at a molecular level, 
thus meaning that the holistic point of view has been somewhat 
overlooked [6]. Indeed, the role of lymphoid organs as the targets 
of infection, the blood and lymphatic circulation, or the intrinsic 
anatomy of the lung (the key organ in this infection) have been 
little studied or ignored. 

SChackerian [7] and Wolf [8] arrived at the same conclusion, 
namely that the most important event for achieving early 
control is fast infection of the hilar lymph node (as production 
of specific lymphocytes is induced earlier), using two different 
experimental models of TB (in mice and macaques). In contrast, 
Bru and Cardona [9] arrived at a slightly different conclusion 
using virtual models, namely that although the time at which 
the immune response is triggered is important, the local 
production of chemokines which can attract those lymphocytes 
to activate the infected macrophages is equally important. This 



supported the idea of granuloma formation being necessary as a 
tool for controlling the infection by making the infectious focus 
"visible" to the specific lymphocytes [10], in contrast to the 
proposal that the granuloma could be a "foe" by allowing the 
bacilli to persist [11]. Observations from other researchers [Ian 
Orme, personal communication] and our group [9,12] also 
support the beneficial role of the granuloma. Thus, no lesions 
can be seen in the lung during the first two weeks of a low dose 
aerosol infection despite the global bacillary load (BL) increasing 
from 10 2 to 10 5 , thereby indirectly demonstrating a "hidden" 
initial phase of the infection [6] . This unicellular infection phase 
appears to last about 10-15 days [9,13] and coincides with the 
time need to detect the first granulomas [7,8,14]. During this 
period, M. tuberculosis grows inside individual infected macro- 
phages that cannot be detected by specific lymphocytes, 
probably because they are not able to produce sufficient 
chemokines to attract them. The existence of this phase could 
also support previous and recent results claiming the detection 
of bacilli in healthy tissue [15,16]. 

In order to study the temporo-spatial relationship between the 
three organs (lung, spleen and lymph node) involved during initial 
M. tuberculosis infection, we designed an experimental study in 
which mice were infected via aerosol (AER) or endovenously (EV) 
and consequently followed-up, subsequendy developing a math- 
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ematical approach to reproduce the results obtained experimen- 
tally and to extract conclusions. 

Results 

Higher Bacillary Doses Trigger an Earlier IFN-y Response 

For consensus, and taking into account the experimental 
results, we considered 20 SFU/250,000 cells to be the threshold 
for onset of the specific IFN-y response. This value was chosen 
because, whenever a significant change was detected in the 
ELISPOT results, it always meant that this value had been 
exceeded. In the EV model, the immune response was first 
detected in the spleen in all cases (Table 1). A comparison 
between the specific responses shows how the immunity mainly 
targets ESAT-6. At the highest dose (10 5 CFU/ml), specific 
lymphocytes are detected first in the lung and then in the 
lymph node, thus reflecting the accumulation of circulating 
IFN-y producing lymphocytes from the spleen before immune 
specific IFN-y response is triggered in the lymph nodes. This 
suggests that some degree of BL is needed for the induction of 
immunity. In the EV model, higher challenge doses lead to 
earlier IFN-y responses. 

In contrast, there is no difference as regards onset of the IFN-y 
response among the organs in the aerosol model: it appears at the 
same time in all of them. 

A Specific Bacillary Load Must Be Reached to Trigger the 
IFN-y Response 

Considering the minimum time required to trigger the IFN-y 
response obtained from the highest challenge model (10 J ), we 
estimated the delay between this trigger and onset of the IFN-y 
response to be 2.5 days. In light of this value, we evaluated the BL 
threshold required to trigger the IFN-y response by analysing the 
level thereof 2.5 days prior to onset of the IFN-y response in each 
experimental model. Specifically, in each case we determined 
which organ had the highest BL 2.5 days before the onset (i.e., 
which the triggering organ was) and which this threshold value 
was. The results showed that the BL level required is relatively 
constant, at around 3.5 loglO (Table 2). This threshold is exceeded 
from the beginning in the highest dose EV model, thus validating 
the hypothesis used in these estimations (i.e., that the 10' dose is 

Table 1. Detection of PPD-specific IFN-y response. 



enough for triggering the IFN-y response immediately after 
infection). 

Control of the BL Is Not Immediate After Onset of the 
IFN-y Response 

Figure 1 shows the evolution of the BL in the different organs as 
obtained in the experimental in vivo models (in black). After the EV 
challenge, the IFN-y response was estimated to appear at days 2.5 
(initial inoculum = 10 5 ), 6.47 (10 4 ) and 1 1.77 (10 2 ) (Table 2); BL 
control occurs subsequently. The sooner the IFN-y response is 
triggered, the sooner BL control occurs, at least in spleen and LN. 
This effect is not seen in lungs. Table 3 shows the delay between 
the appearance of the IFN-y response and the effective control of 
the BL, per organ. 

The BL in an organ does not appear to be crucial as regards the 
control of bacillary growth, as this occurs over a wide range from 
10 5 to 10 7 CFU/ml depending on the initial infective dose. 

In the aerosol-infected model, BL control in the lung, which is 
the organ with the highest bulk, appears at the same time as the 
immune response. However, and as is the case for the EV model, 
BL control in the spleen is achieved 5 days after the onset of the 
IFN-y response (i.e. at day 16.43). This value confirms that, 
although a low BL in the spleen is able to trigger and attract the 
IFN-y response, it is not capable of effectively stopping bacillary 
growth. 

The Mathematical Model Designed Reproduces the 
Experimental Data and Highlights the Importance of 
Considering the Spatial Factor 

The in silico model built provided results which agreed well with 
experimental measurements in terms of time and levels (Figure 1). 
10 4 and 10 5 series were used for parameter estimation (i.e., for 
fitting the model to experimental data), and 10 2 and AER series 
where used for validation of the model. However, BL control was 
achieved earlier in the mathematical model than in the in vivo 
models under all conditions. 

Discussion 

The aim of this study was to evaluate the bacillary load and 
specific IFN-y responses in lung, spleen and lymph node, after 
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Day of detection of the PPD- and ESAT-6-specific IFN-y responses in the three different organs obtained for all experiments performed. 
doi:1 0.1 371 /journal.pone.01 00830.t001 
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Table 2. Onset of the specific IFN-y response. 





APPEARANCE OF THE SPECIFIC IFN-y RESPONSE 


EV exp (depending on the Inoculum, in CFU) 




AER exp 


10 2 10* 


10 5 




Organ of appearance of the spleen spleen 
specific IFN-y response (>20 SFU) 
and day of appearance 
of such response 


spleen 


spleen and lymph node 



11.7±0.4days 6.4±0.9 days 2.5±0.5 days 16.4±2.3 days 



Organ that triggers the spleen spleen spleen lymph node 

specific IFN-y response 

and its BL 2.5 days before 

appearing the specific IFN-y 

response 

3.41 ±0.21 log CFU 3.77±0.17 log CFU 4.66±0.14 log CFU <3.61 ±0.23 log CFU 



First row: estimated onset of the specific IFN-y response taking into account the organ and day on which it appears. The moment at which the 20 SFU threshold would 
have been surpassed is linearly interpolated between two consecutive experimental points (last point with SFU<20 and first point with SFU>20). Second row: organ 
and BL of the organ that would have triggered the IFN-y response 2.5 days before its estimated appearance (i.e., organ with the highest BL at that moment}. 
doi:1 0.1 371 /journal.pone.01 00830.t002 



aerosol or intravenous infection; and to develop a mathematical 
approach with the data obtained. 

To date, a paramount role has been assigned to rapid infection 
of the lung lymph node as regards onset of the immune response 
and better infection control in the lung [7,8]. The results showed 
that a critical BL has to be reached to trigger the immune 
responses (at least in terms of IFN-y responses) and that as highest 
is the BL the earlier this response is triggered, and yet 
(paradoxically) it is not sufficient to stop immediately the BL 
growth. 

The dual role of lymph nodes and spleen as providers of 
effector immune cells and as infectious foci intrigued us, and 
anatomy is one of the possible answers. An infected macrophage 
or dendritic cell could reach the trabecular sinuses, thereby 
initiating infection of the lymph node in an anatomically 
different location in the cortex or the medullar sinus, where 
antigens are presented to B and T cells, respectively [17]. If this 
is the case, local infection in the lymphoid tissues would not be 
directly controlled by the local IFN-y response triggered, which 
will be spilled to the efferent capillary lymph vessels and 
therefore to pulmonary circulation thus to the right heart and 
lungs [18], where it is attracted and further captured by 
pulmonary infectious sites. Likewise, the spleen, which special- 
izes in retaining and destroying old erythrocytes [19], has two 
compartments: the red pulp, where an infection is more likely to 
occur, and the white pulp, in which the antigens are presented 
and specific lymphocyte proliferation induced [19]. 

SOur initial hypothesis was that this was possible because 
lymphatic organs could host an infected site while still being 
able to develop a specific IFN-y response against this infection 
(although not benefiting from it) for anatomical reasons. Indeed, 
the results of our mathematical model support this. Because of 
the anatomy of lymphatic tissues, the lymphocytes would follow 
the lymphatic and blood circulation and be further attracted by 
other infectious foci. In the case of a low dose aerosol at day 2 1 
post-infection, when the bacillary load in the lung is very high, 
a lymphocyte generated in the spleen will reach the cava vein 
and then progress further to the right atrium and right 
ventricle, eventually reaching the lung. At this point, because 
of the high infection, the lesions are likely to attract this 



lymphocyte. If this does not occur, it will continue into the left 
atrium and ventricle to circulate throughout the rest of the 
body, perhaps being attracted by the infected lymph nodes or 
even reaching the spleen. This dynamic was validated by the 
mathematical model. 

Interestingly enough, even in this situation it is not clear 
whether these lymphocytes will be able to activate the infected 
macrophages, as a pool of attraction (i.e. a local focus of 
chemokine production) is needed, suggesting the granuloma 
allows the attraction of specific lymphocytes [10]. Considering 
that the initial infection usually takes place inside the alveolar 
macrophage (in the alveolar space), where immunological 
surveillance is almost impossible, favoring reinfection, the lack 
of pool of attraction is a problem in terms of infection control 
in this timeframe [9,20]. 

We are aware that the host immune response against 
tuberculosis involves a complex combination of cytokines, 
chemokines and signaling molecules other than the IFN-y 
responses, and our study was only focused on PPD and ESAT-6 
specific IFN-y responses. However, as these have been widely used 
as a marker of tuberculosis infection even if its limitations are 
commonly known and accepted, we did considered that they could 
fairly be representative of the host immune response for the study 
purposes. Moreover, the number of animals involved was also 
relevant for this decision, as to add more essays to explore other 
would have meant to increase this number in detrimental of ethics. 
On the other hand, we tried to keep the parameters and premises 
of the in silico model at minimum, in order to keep it as simple as 
possible. 

Another limitation of our study is that it is based on an 
experimental mice model, and mice are far from humans. 
However, in terms of anatomy of the organs studied and dynamics 
of blood circulation, the results obtained in mice can be totally 
extrapolated to humans, at least those related with the qualitative 
dynamics. The quantitative results (for instance, the delay between 
trigger and onset of IFN-y response, or the BL threshold for such 
triggering) may slightly change from mice to human, and thus the 
mathematical model should be specifically parameterized for the 
second case. 
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Figure 1. Output of the in silico model and the corresponding experimental results. The output of the in silico model (white symbols) are 
represented considering EV infections of 10 2 , 10 4 and 10 5 CFU/mL and an AER infection. The experimental results are represented by black symbols. 
The results are shown individually for each organ (lungs (circles), spleen (squares) and lymph node (triangles)). Ten repetitions of each simulation 
were performed by changing the random seed. Each experimental point corresponds to an individual mouse. 
doi:1 0.1 371 /journal.pone.01 00830.g001 



The mathematical model designed to simulate the complex 
anatomical relationships between the different infected organs 
reproduces the general dynamics of the experimental results. As 
in the generation of any other in silico model, assumptions have 
to been made and the parameters chosen, something we did 
trying to keep it simple and in accord with experimental 
observations. The only point where it differs from them is on 
the delay on controlling the BL, as it only takes into account 
the global BL of each one without considering the local 
organization of the infected cells, and cannot incorporate the 
spatial dimension which might be important for attracting 
specific lymphocytes to each infectious focus. This may be the 
reason why faster control of the BL is achieved after onset of 



the IFN-y response than in the in vivo model, therefore showing 
the importance of this factor. 

In conclusion, although lymphoid organs appear to be 
organized as a two compartment structures, thus favoring infection 
control in non-lymphoid organs, the the local attraction of specific 
lymphocytes to the infected sites is crucial as regards controlling 
the BL in these organs and on the infection control. 

Overall, these data suggest that the key to subversion of the 
immune response by M. tuberculosis is not related to the ability of 
the bacilli to slow its induction. Indeed, our findings show that the 
speed of the immune response is not related to the degree of 
infection control as the ability of infected cells to attract specific 
lymphocytes appears to be the most critical factor as regards 
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Table 3. Delays between the appearance of the IFN-y response and the effective control of the BL, per challenge dose and route, 
and per organ. 







Spleen 


Lymph node 


Lungs 


EV exp 10 2 


Day of estimated appearence of the IFN-y response (1) 


11,7 






Day of detected control of BL 


17 


17 


14 


Delay between appearence of the IFN-y response and BL control (days) 


5,3 


5,3 


2,3 


EV exp 10 4 


Day of estimated appearence of the IFN-y response (1) 


6,4 






Day of detected control of BL 


14 


14 


11 


Delay between appearence of the IFN-y response and BL control (days) 


7,6 


7,6 


4,6 


EV exp 10 s 


Day of estimated appearence of the IFN-y response (1) 


2,5 






Day of detected control of BL 


11 


11 


14 


Delay between appearence of the IFN-y response and BL control (days) 


8,5 


8,5 


11,5 


AER exp 


Day of estimated appearence of the IFN-y response !1) 


16,4 




16,4 


Day of detected control of BL 


21 


21 


16 


Delay between appearence of the IFN-y response and BL control (days) 


4,6 


4,6 


=>0 



Deviations shown in Table 2 

These delays are evaluated by subtraction between the day when BL control is experimentally detected at each organ and the day of the estimated appearance of IFN-y 
response. 

doi:1 0.1 371 /journal.pone.01 00830.t003 



benefitting from the immune response, at least in terms of IFN-y 
responses. 

This might have a conceptual impact when it comes to 
designing new prophylactic strategies, especially new vaccine 
candidates as, if they are able to trigger a rapid IFN-y response 
(or even -if we hypothesize- a rapid Th 1 immune response), will 
this be sufficient to cope with the infection? The answer to this 
is probably not, as this enhanced IFN-y response will not know 
where to act as no infectious foci will be generating sufficiendy 
strong attractive signals. In light of the results of our study, we 
therefore suggest that other strategies focusing on therapeutic 
vaccines or those generating humoral responses, for example, 
should perhaps be explored. 

Materials and Methods 

Experimental Data 

Mice. 6-8-Week-old female C57BL/6 specific-pathogen-free 
(spjj mice were obtained from Harlan Laboratories (Sant Feliu de 
Codines, Catalonia, Spain). The animals were shipped under 
suitable conditions, with the corresponding certificate of health 
and origin. Upon arrival, the mice were kept under controlled 
conditions in a P3 high security facility with sterile food and water 
"ad libitum". 

Ethics. All animal procedures were approved and supervised 
by the Animal Care Committee of the Germans Trias i Pujol 
University Hospital and by the Department of Environment of the 
Catalan Government. Mice were weighed and checked every week 
following a protocol that monitored weight loss, apparent good 
health (bristled hair and wounded skin) and behaviour (signs of 
aggressiveness or isolation). Mice were euthanized with isoflurane 
(inhalation excess) in all cases in order to avoid any suffering. 

Experimental plan. Two routes of infection were evaluated: 
the aerosol and the intravenous route. A total of six animals per 



experimental group were included in each experiment. The AER 
experiment was run twice, whereas the EV experiment was run 
only once due to the already high number of mice involved for the 
three different inoculation doses being evaluated (10 2 , 10 4 or 10 5 
Colony Forming Units (CFU) per mL). 

The animals were sacrificed at established sequential time- 
points. Different timepoints were established for AER and the EV 
exp based on our previous experience with experimental TB 
infection. In the aerosol experiment, animals were sacrificed on 
each day from day 10 to 22 post-infection, whereas in the EV 
experiment mice were sacrificed on days 0 and 2, 5, 8, 1 1, 14, 17, 
20, 28 and 42 post-infection. 

Infection. M. tuberculosis strain H37Rv Pasteur was grown in 
Proskauer Beck medium. In the AER experiment, animals were 
placed in the exposure chamber of an airborne infection apparatus 
(Glas-col Inc., Terre Haute, IN, USA) for infection with a low dose 
of M. tuberculosis. Nebulization provided an approximate uptake of 
20-50 bacilli by mice lungs. 

In the EV experiment, a total volume of 0.2 mL containing 10 2 , 
10 4 or 10 5 CFU per mL (depending on the experimental group) 
was inoculated into the tail vein of each mouse. 

Evaluation of the infection. The bacillary load (BL) and 
specific IFN-y responses measured were the parameters chosen to 
evaluate the infection. Both parameters were measured in lung, 
spleen and lymph nodes, which were extracted and mechanically 
disrupted in order to obtain tissue homogenates. Lymph nodes 
samples could occasionally not be harvested at the earliest 
timepoints as they were too small to be readily recognised, 
therefore no results were recorded in such cases. 

The BL was determined by culturing the samples on 
Middlebrook 7H11 agar plates (Bennex Ltd, Shannon, Ireland) 
at 37°C for 21 days. Visible colony forming units (CFU) were 
counted four weeks later, with data being recorded as the log of the 
total number of CFUs recovered per organ. 
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Specific cellular IFN-y responses against the antigens ESAT-6 
(Lionex Diagnostics & Therapeutics, Braunschweig, Germany), 
which was used at a final concentration of 5 \lg/ mL, and PPD 
(Statens Serum Institute, Kobenhavn, Denmark), which was used 
at a final concentration of 10 |0.g/mL, were measured using an 
IFN-y ELISPOT kit (BD Bioscience, San Diego, CA) after plating 
250,000 cells from the tissue homogenates. The ELISPOT assay 
was further developed and read according to the manufacturer's 
instructions. The results were expressed as IFN-y secreting cells as 
Spot Forming Units (SFU) per 250,000 cells. 

The In Silico Model 

An in silico model based on several premises was designed in order 
to test our hypothesis regarding the in woo dynamics of the infection. 
It was fitted to and validated with experimental results and 
subsequendy further used to extract some conclusions about how 
the flows between the different organs influence the infection. In 
order to build this model, different variables, processes and 
parameters were taken into account. These are described below (SI). 

Model Variables 

Those quantities that vary with time and represent the dynamics 
of the system were termed variables. According to the experi- 
mental measurements, these variables are the BL and the level of 
IFN-y response expressed as IR. 

BLi(t), BL s (t) and BLi n (t) are the bacillary loads in lung, spleen and 
lymph nodes at time t, respectively, and IRi(t), IR,(t) and IRi n (t) the 
immune responses in lung, spleen and lymph nodes at the same time 
t, respectively. The dynamics of these variables are governed by the 
model processes described below. The output of the model should 
agree with the experimental behavior of these variables. 

Model Processes 

The model processes are those that can make the variables 
change, and are based on assumptions decided by consensus after 
checking different possibilities and their effects on the model 
output, that should be in accord with experimental data. Three 
such processes were considered: local bacillary growth, the trigger 
for the specific immune response and the flows between the three 
different organs studied. 

Local bacillary growth. In order to model the bacillary 
growth in each organ (lung, spleen and lymph nodes), the effective 
net growth arising from replication inside macrophages and 
bacterial death was taken into account. This growth, which is 
represented by the classic logistic model, is modulated by the IFN- 
y response activity (a IR ) as follows (Eq. 1). 

dBLi 

— ; — =u-max 

(1) 

(0,BL,(t) - am-IRM) ' ( 1 - BL^v) ' = Z " ? ' ln 



organ. This is thought to occur after the lysis of infected 
macrophages, when bacilli are released into the extracellular 
environment. As such, migration of bacilli from one organ to 
another is linked to the discrete intervals at which macrophages 
burst. Consequendy, this process is modelled as discrete migrations 
at certain time intervals corresponding to multiples of T lysis after 
initial infection of the corresponding organ rather than as a 
constant flow. The parameter Ti ysis represents the mean time that 
a macrophage takes to burst once infected. 

The physiological flows (determined by the blood circulation 
dynamics) important (i.e., numerically relevant to the results) for 
the dynamics of the model were chosen according to the 
principle of Occam's razor [21]. That is, all such flows were 
initially incorporated into the model but only those that were 
numerically relevant were kept. In the EV case, initial 
concentrations in the spleen and lungs were similar. If lungs 
were not initially infected (low concentration), infection came 
from the spleen. Some of the bacilli coming from the spleen 
could not be detained in the lungs and arrived directly at the 
lymph nodes without being previously phagocytosed by an 
alveolar macrophage. In the AER case, the main flow was from 
lungs to lymph nodes and then on to the spleen. Table 4 shows 
the flows between organs considered. 

Macrophage lysis can be considered to be a synchronic 
phenomenon at the beginning of the infection (i.e., infected 
macrophages burst approximately simultaneously), although it 
becomes desynchronized as time goes by. As such, bacillary 
migration should be considered to be discrete at the beginning 
of the infection and as continuous at some point thereafter. 
However, this phenomenon is only significant at the beginning 
of the infection as, when the bacillary load is high, the 
migratory effect between organs becomes irrelevant. 

Trigger and dynamics of the specific IFN-y 
response. The experimental results for the EV experiment show 
that the IR depends on the inoculated dose and, therefore, that it 
should also depend on the bacillary load: while BL remains small, 
no IR would be activated. Once a certain threshold is surpassed, an 
IR would be triggered. Another important experimental observa- 
tion is that although very high inoculated doses (10 5 CFU/mL) 
should be sufficient to activate the IR, IR measurements take some 
time to become positive. This suggests that, once the IR is triggered, 
there is some delay until it appears. These two observations are 
incorporated into the model using the parameters BL, (lrt , s and T delay , 
which represent the BL threshold to trigger the IR and the time 
delay for the detection of this response, respectively. 

Once the IR is activated and generated in a lymphatic organ, 
lymphocytes are thought to enter the blood stream in order to 
reach the infected organs. As such, the model should consider the 
specific IR level in blood (IRbiood) as an auxiliary variable. This 
variable increases with lymphocyte production at the rate y, driven 
by BL thres and T delay , and decreases due to lymphocyte lysis (with a 
rate CO) or when the lymphocytes are retained by an infected organ, 
as shown by Eq (2). 



Thus, Eq. (1) shows the time (t) evolution of BL at each organ (1 dIRbi oo d(0 
for lungs, s for spleen and ln for lymph nodes), as modulated by the ^ 

IFN-yresponse level (IR) in the corresponding organ (i = l, s, ln). dIR(t) 

The parameters involved are the bacillary effective growth rate (u), =H i (BL i (t—T dL ,i ay ))-y i — IR h i ood (tya)— ^ — -j— 

the maximum bacillary load per unit volume (BL max ), the effective j=l,s,in ' 

organ volumes (V;, i = l, s, ln), and the IFN-y response activity , ■. 

(am). The latter represents the number of bacilli that are wUh H . (x) = \ * X< ,hres I and i = sM 

inactivated per unit of specific IFN-y response. I I if x> BL t i, re s J 

Bacillary flows between organs. The bacilli grow locally in 
each organ but can be released into circulation and reach another 
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Table 4. In silico model of bacillary migration between organs, where BLj(t) is the bacillary load in spleen (i = s), lungs (i = l) and 
lymph nodes (i= In) at time t. 



Bacilli flow 



Mathematical expression 



Spleen (s) to 
lungs (I) 

Spleen to 
lymph nodes (In) 

Lungs to 
lymph nodes 

Lymph 

nodes to spleen 



BL, (t)*-BLj (i) + ii s ^if(t-t il ,fJ)modTi y „ = 0 
BLi„{t)^BLi„{t) + C s ^ h y{t- ti„ fj ) modTi ym = 0 
BLi„(l)^BLi„(t) + l, -,,„ if (t - ttf, i ) modT h ,„ = 0 
BL,{t)^BL,{t) + Ci a ^Jf(t-ti,, f , ln )modTi, les = 0 



The auxiliary variables t inf s , t in fj and t infl „ represent the time that spleen, lungs and lymph nodes have been infected. The parameters involved are detailed in Table SI. 
doi:1 0.1 371 /journal.pone.01 00830.t004 



The retention of an IR by an infected organ is assumed to 
depend on the bacilli growing in that organ (i.e., it does not 
depend on the bacillary load of the organ, BL ; (t), but on BL ; (t)- 
aiR-IR;(t), as shown in Eq (1)). Other parameters involved in 
lymphocyte transport include the blood flow F blood (volume of 
blood that enters the left ventricle of mouse heart per unit time) 
and the fraction of this flow that enters each of the organs 
considered (Qj with i = 1, s, In). Eq (3) shows the equation for the 
model governing lymphocyte transport to these organs. 



dlR; 



= IB-blood (*)' Fblood' Qr 

Max(0,BLj(t) - a IR -IRi(t)) 



(3) 



--Is An 



The model presented above for the dynamics of the IR is 
extremely simple but attempts to represent the overall dynamics of 
IR by taking into account the most relevant processes. 

Parameter Estimation 

Parameters are those quantities that determine the result of the 
processes studied, for instance, the effective bacillary growth rate 
(H). Their values must be directiy extracted or estimated from 
experimental data and the literature or, in the worst case, inferred 
by fitting the model output to experimental results. Table 4 shows 
the list of parameters of the model, as well as their values and how 
these values have been fixed. 

The duplication time was set to 24 hours, which is in 
accordance with the literature [22] and corresponds to a specific 
growth rate of 0.029 h _1 . This value is compatible with the 
absolute maximum apparent growth rate observed experimen- 
tally, in particular in the EV experiment with low inoculum (10 
CFU/ mL) in spleen. In these measurements, if we consider the 
interval between days 2 and 10 (Fig. 1) we obtain a duplication 
time of 1.39 days (R 2 = 0.999). This value apparendy shows a 
slower growth but is distorted by local factors such as bacillary 
flows. The parameters T de i ay and BL thres were inferred from 
experimental measurements, as detailed in Section 2.2. The 
extinction rate of specific lymphocytes (oo) was estimated using 
the expression co= ln2/t 1/2 , where t 1/2 is the time at which the 
specific IFN-y response is reduced to half of its maximum value. 
Reference [9] provides an order of magnitude of 3 hours for 
this value, which is in accordance with the half-life of specific 
lymphocytes in the absence of antigens. The percentage of 



blood that enters each organ, Qj, was chosen on the basis of the 
irrigation for each organ. 

Finally, the remaining parameters were estimated by fitting the 
model simulation results to EV experimental results for inocula of 
10* and 10 J CFU/mL, only. To do so, an objective function was 
defined (Eq. 4): 



f(BL)= 

\<k<N 



BL mo d(t/ c ) — BL exp ( tk ) 



BL^Jt/c) 



\ 



lung 



BL mo( a (tk) — BL cxp ( tk ) 



BL exv (tk) 
BL mo( j(tk) — BL cxp (tk) 



BL Exp (t k ) 



spleen 



lymph 
node 



(4) 



The model was implemented in the computing language C++ in 
order to solve those equations that had no analytical solution, 
which were solved using Euler's numerical method (controlling the 
instabilities). An exhaustive search of the set of parameters that 
minimized this objective function was performed in order to 
minimize differences between model output (mod) and experi- 
mental measurements (exp). To do so, two series of simulations of 
the mathematical model were performed, one for each initial 
concentration. Each series of simulations consisted of several 
repetitions with the same initial conditions, changing the random 
seed in order to reproduce experimental variability. 

The initial bacillary load values were taken from experimental 
measurements by extrapolating data from day 2. These values 
were BL { (0) =10 CFU, BL S (0) =270 CFU and BL ln (0) = 1 CFU 
for the former and BL, (0) = 2500 CFU, BL S (0) = 1 1000 CFU and 
BL ln (0) =10 CFU for the later. The outputs obtained from the 
model are shown in Figure 1. As can be observed, the initial 
increase of infection in several organs and the plateau once the 
specific IFN-y response is activated are well reproduced. The best 
values for the parameters estimated are shown in Table 4. 

Model Validation 

A second set of simulations was performed using the parameters 
estimated in the previous series. In this case only the initial 
conditions were modified in an attempt to reproduce AER and 
low concentration EV experiments without varying the model 
parameters. The initial infections were fixed as 10 bacilli in lungs 
(AER) and 5 bacilli in spleen (low concentration EV) with the same 



PLOS ONE | www.plosone.org 



7 



June 2014 | Volume 9 | Issue 6 | e100830 



An Earlier IFN-y Response, Not Sufficient to Control M.tuberculosis 



extrapolation method (see previous section). As shown in Figure 1, 
the model results were in agreement with experimental measure- 
ments. In the AER case, the coincidence between model and 
experimental results is high, although a slight underestimation of 
the plateau values can be observed. The low concentration EV 
simulations reproduce the infection evolution in spleen and the 
high variability of measurements in lungs and lymph nodes, which 
is also in agreement with experimental results. 

The underestimation of the plateaus seen in most of the 
simulations suggests that although the model considers the effect of 
the specific IFN-y response to be immediate, this might not be the 
case in vivo (as already suggested by the experimental results). 



Supporting Information 

Table SI Parameters of the in silico model. The 

Parameters of the in silico model described in Equations (1), (2) 
and (3) and in Table 3 have been included in this table, as well as 
their description, source and values. 
(PDF) 
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